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! Abstract 

' The energy spectra of ultra high energy cosmic rays reported by the AGASA, 

Fly's Eye, Haverah Park, HiRes, and Yakutsk experiments are all shown to be in 
i agreement with each other for energies below 10^'' eV (after small adjustments, 

I within the known uncertainties, of the absolute energy scales). The data from 

HiRes, Fly's Eye, and Yakutsk are consistent with the expected flux suppression 
, above 5 x 10^^ eV due to interactions of cosmic rays with the cosmic microwave 

I background, the Greisen-Zatsepin-Kuzmin (GZK) suppression, and are inconsistent 

CN ' with a smooth extrapolation of the observed cosmic ray energy spectrum to energies 

^ ! > 5 X 10^^ eV. AGASA data show an excess of events above 10^'' eV, compared to 

[ the predicted GZK suppression and to the flux measured by the other experiments. 
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X 



1 Introduction 



We analyze the observed spectrum of ultra-high energy cosmic rays. We find 
two main results: (i) The energy spectra reported by the AGASA, Fly's Eye, 
Haverah Park, HiRes and Yakutsk experiments are all in good agreement for 
energies below 10^0 eV, and (ii) All the data are consistent with a GZK sup- 
pression except for the AGASA points above 10^° eV. Our principal conclusion 
from these two results is that standard physics, including the GZK suppres- 
sion, is sufficient to explain all of the existing data on UHE cosmic rays. 

For any theoretical model in which the GZK suppression is present, the as- 
sumed intrinsic spectrum produced by the UHE cosmic-ray sources influences 
the energy spectrum predicted by the model. Our conclusion that the data 
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Fig. 1. Currently available data on the highest energy cosmic rays. The quantity J is 
the differential energy flux, dcp/dE, per unit time per unit area per steradian. Panel 
(a) shows the data as published by the five experimental collaborations: AGASA 
[1], Fly's Eye [2], Haverah Park [3], HiRes [4], and Yakutsk [5]. Panel (b) shows the 
data after adjusting the absolute energy calibrations of the various experiments so 
as to bring the results from the different experiments into agreement at 10^^ eV. 
For specificity, the Fly's Eye absolute energy scale was adopted as the standard. 
The fractional shifts in absolute energy scale, AE/E, shown in the figure, are all 
well within the published systematic errors in the energy scale. 

are consistent with a GZK suppression implies that the observed spectrum is 
consistent with model predictions for a plausible intrinsic energy spectrum. 
In particular, we show that the observed spectrum is consistent with that 
expected for a GZK suppression of the flux produced by a simple cosmologi- 
cal distribution of sources, each source producing high energy protons with a 
spectrum dN/ dEp oc E'"^ characteristic for collisionless shock acceleration. 

Before entering into any details, we will summarize and compare in this intro- 
duction the data that are available from different collaborations that measure 
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the spectrum of ultra high energy cosmic rays. 

1.1 Summary of available data 

Figure 1 is a " Before- After" figure of the currently available data on the high- 
est energy cosmic rays (energies > 10^^ eV). In Figure la (the "Before" ver- 
sion of the figure), the data are plotted, together with their flux error bars, as 
they have been pubhshed by the flve experimental collaborations: AGASA [1], 
Fly's Eye [2], Haverah Park [3], HiRes [4], and Yakutsk [5]. The Haverah Park 
data have recently been re-analyzed using modern numerical simulations of 
air-shower development [3]. The reanalysis resulted in significant changes of 
inferred cosmic-ray energies compared to previously published results ([6] and 
references quoted therein). The data points for the Haverah Park measure- 
ments that are shown in Fig. 1 are based on this improved analysis, which is 
available only at energies < 10^^ eV^ . 

The most striking feature of Figure la is that the experimental results differ 
greatly among themselves (by factors ~ 2) even in the region 10"^^ eV < E < 
2 X 10^^ eV, where the quoted error bars from each experiment are very small. 
In addition, the higher AGASA flux reported above 2 x lO^'' eV stands out 
above the scatter in the different experimental measurements. 

Figure lb (the "After" version of our " Before- After" flgure) shows a dramat- 
ically different representation of the available data. With small adjustments 
in the absolute energy scales, all of the measured fluxes are seen to be in 
agreement at energies below 10^° eV. In constructing Figure lb, we have ad- 
justed the absolute energy calibrations within the error bars published by the 
experimental collaborations. We chose the shifts so as to bring the different 
measured fluxes into agreement at 10^^ eV. The energy shifts can be accom- 
plished in five equivalent ways, depending upon which one of the five energy 
scales is unaltered. For Figure lb, the Fly's Eye energy scale was unaltered 
and we adjusted the AGASA energy scale by -11%, Haverah Park by +15%, 
HiRes by +7.5%, and Yakutsk by -19%. All shifts are well within the published 
systematic errors. 

Figure 1 illustrates visually our two main points. First, all of the currently 

^ A single flux point at '-^ 7 x 10^^ cV is shown in Figure la, but this point is 
based on a preliminary analysis of 4 events that are chosen by different cuts than 
those applied for the lower energy data. The energy uncertainty for the point at 
~ 7 X 10^^ cV is significantly larger than the estimated uncertainties for lower 
energy points [3]. Therefore, the point at '-^ 7 x 10"*^^ eV is shown in Figure la only 
for completeness; it is not used elsewhere in our analysis because the Haverah Park 
collaboration has described this point as preliminary. 
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available data on high energy cosmic rays are in agreement within their quoted 
errors for energies between 2 x 10^^ eV and 10^'' eV. Second, three of the four 
data sets available above 10^^ eV, HiRes, Fly's Eye, and Yakutsk, all show 
evidence for a turnover of the energy spectrum for energies above 5x10^^ eV. 
This turnover, we shall show later, is highly significant statistically and is 
consistent with what one would expect from a simple model that includes the 
GZK effect. Above 10^° eV, the reported AGASA fluxes are higher than the 
fluxes measured in other experiments. It is these high AGASA fluxes alone 
that have led to the widespread impression that measurements of ultra-high 
energy (UHE) cosmic rays (energies > 10^^ eV) do not show evidence for a 
GZK effect. 

1.2 What does it all mean? 

What can one make of the results shown in the Before-After Figure 1? There 
are two simple possibilities. First, the excellent agreement shown in Figure lb 
among the different experiments could be accidental. According to this inter- 
pretation, the small adjustments made in the energy scales are not physically 
motivated and the real situation is somehow much more complicated. It is just 
a fluke that all of the adjusted energy spectra line up together so well below 
10'^'' eV. This interpretation is certainly possible. In the present paper, how- 
ever, we shall choose a different interpretation of Figure lb. We shall suppose 
that the excellent agreement of the adjusted energy spectra reveals a good 
approximation to the true shape of the UHE cosmic ray energy spectra. We 
shall now explore the consequences of this assumption. 

We stress that the distinction between the two possibilities for interpreting 
Figure lb can only be settled by a new generation of precise and high statis- 
tics measurements of the UHE cosmic ray spectrum. Fortunately, the Auger 
experiment, currently under construction [7], is expected to provide the nec- 
essary precision and statistics. The Telescope Array experiment [8], currently 
under planning, may also provide similar precision and statistics. 

We first describe in Section 2 the model we use and then in Section 3 we com- 
pare the model predictions with observations of the UHE cosmic ray energy 
spectrum. We summarize our main conclusions in Section 4. 



2 A simple two-component model 

We describe in this section a simple two-component model for the energy 
spectrum of the highest energy cosmic rays. The Galactic component is taken 
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from observations of the Fly's Eye group. The two input parameters for the 
extra-galactic component (the rate of energy deposition in cosmic rays and 
the shape of the initial spectrum) were originally suggested by the idea [9-11] 
that gamma-ray bursts (GRBs) are the source of UHE cosmic rays. However, 
any postulated cosmologically distributed source of cosmic rays with a similar 
energy production rate and energy spectrum (Eq. 1 and 2 below) would yield 
agreement with the observations. 

In order to avoid the risk of being misled by "curve fitting", we use the same 
theoretical model that was discussed in 1995 [12]. We assume that extra- 
galactic protons in the energy range of 10^^ eV to 10^^ eV are produced by 
cosmologically-distributed sources at a rate 

de 

— fti 3 X lO^Vg Mpc-3 yr-\ (1) 

with a power law differential energy spectrum 

^ocV. "-2. (2) 



Wc shall refer to this energy spectrum as "the extra-galactic component" in 
order to emphasize that the fit to the data is generic, independent of the type of 
source that generates the assumed energy and spectrum. An energy spectrum 
similar to the assumed energy spectrum, Eq. (2), has been observed for non- 
relativistic shocks [13] and for relativistic shocks [14] shocks. This power law 
is produced by Fermi acceleration in collisionlcss shocks [13], although a first 
principles understanding of the process is not yet available (see, e.g. Ref. [15] 
for a discussion of alternative shock acceleration processes). 

We can use Eq. (1) and Eq. (2) to obtain a value for the cosmological rate, 
E^dN/dEp, at which energy in the form of high energy protons is being pro- 
duced. Integrating EpdN/dEp between 10^^ cV and 10^^ eV and setting the 
result equal to the value given in Eq. (1), wc find the proportionality constant 
in Eq. (2). Thus E^dN/dEp ^ 0.7 x lO^^gMpc'^ yr'^ 

Energy losses due to pion or pair production are included in the transport 
calculations in the usual continuous approximation. Energy loss due to the 
cosmological redshift is significant for energies < 5 x 10^^ eV, and is also 
taken into account. The choice of cosmological model is unimportant for cos- 
mic ray energies above 10^^ eV, which is the region of interest. We assume, 
for definiteness, a flat universe with Qm = 0.3 and = 0.7, and Hubble 
constant Hq = 65km/sMpc. For consistency with our earlier derivation of 
the upper bound on neutrino fluxes that follows from the observed cosmic 
ray spectrum [16], we assume that the source density evolves with redshift z 
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like the luminosity density evolution of QSOs [17], which may be described as 
f(^z) = (1 + z)" with a ^ 3 [18] at low redshift, z < 1.9, f{z) = Const, for 
1.9 < z < 2.7, and an exponential decay at z > 2.7 [19]. This functional form 
of f{z) is also similar [17] to that describing the evolution of star formation 
rate [20] , and also beheved to describe the redshift evolution of GRB rate (see, 
e.g. [21] for review). As mentioned above, the choice of redshift evolution does 
not affect the spectrum above 10^^ eV. 

The cosmic-ray spectrum flattens at ~ 10^^ eV [2,1]. There are indications that 
the spectral change is correlated with a change in composition, from heavy to 
light nuclei [2,22,23]. These characteristics, which are supported by analysis of 
Fly's Eye, AGASA and HiRes-MIA data, and for which some evidence existed 
in previous experiments [6], suggest that the cosmic ray flux is dominated at 
energies < 10^^ eV by a Galactic component of heavy nuclei, and at UHE by 
an extra-Galactic source of protons. Also, both the AGASA and Fly's Eye 
experiments report an enhancement of the cosmic-ray flux near the Galactic 
disk at energies < 10^^'^ eV, but not at higher energies [24]. 

We therefore add an observed Galactic component. 



to the extra-galactic spectrum component given in Eq. (2). The shape of the 
energy spectrum of the Galactic component, Eq. (3), was derived by the Fly's 
Eye collaboration [2] . 

The observed E~'^-^ spectrum between 1 x 10^^ eV to 5 x 10^^ eV is, in this 
model, the combination of two different source spectra. First, the cosmological 
distribution of sources generates an E^"^ spectrum (see Eq. 2), which energy 
losses due to interactions with the CMB steepen to an observed spectrum that 
is a bit shallower than Second, the Fly's Eye fit to the Galactic heavy 

nuclei component makes a small contribution at energies > 1 x 10^^ eV and 
is steeper than E"^-^. We will now compare the model spectrum produced by 
these two sources with the cosmic ray observations. 



3 Compcirison of Model with Cosmic Ray data. 



Figure 2 compares the model prediction with the data from the AGASA [1], 
Fly's Eye [2], Hires [4], and Yakutsk [5] cosmic ray experiments. In order to 
demonstrate that our results are insensitive to the choice of absolute energy 
scale, we present results for three different choices of the absolute energy scale: 
adopting the Fly's Eye, the AGASA or the Yakutsk energy calibration. The 
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Fig. 2. Model versus data. The solid curve shows the energy spectrum derived 
from the t'wo-component model discussed in Section 2. The dashed curve sho'ws the 
extra-Galactic component contribution. The "No GZK" curve is an extrapolation of 
the E'"^'^^ energy spectrum derived for the energy range of 6 x 10^^ eV to 4 x 10^^ eV 
([25]; see text) . Three choices of the absolute energy scale are illustrated. 
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three best fit (solid) curves correspond to energy generation rates (see Eq. 1) 
of de/dt = {2.5, 3.0, 3.5} x lO'^^erg/Mpc^yr and spectral indices, (see Eq. 2) of 
n = {-2.2, -2.1, -2.0} for the {Fly's Eye, AGASA, Yakutsk} energy scales, 
respectively. 

3.1 Good agreement below 10^° eV 

The model predictions are in good agreement with the data of all experiments 
in the energy range 10^^ eV to 10^° eV, a region in which the extra-galactic 
component is predicted to be dominant. Since the Fly's Eye representation of 
the Galactic component is intended to describe the lower energies, it is not 
surprising that the model results are also in good agreement with the observed 
spectrum for energies below 10^^ eV. 

A "X" by-eye" comparison of the model to the data shown in Figures 2 appears 
to indicate a quantitatively good fit, but we cannot simply compute a formal 

fit due to the uncertainties in the absolute energy calibration. Instead, we 
compare the variance of model predictions from the combined AGASA, Fly's 
Eye, HiRes and Yakutsk data sets (^Modei) with the variance of AGASA, Fly's 
Eye, Haverah Park, and Yakutsk data sets from the HiRes data set (s^j^). 
Let SModei = -^"^ T,i,j{nij - Ui j ^ModeiY /riij, Model where nij-,Modei is the predicted 
average number of events in the i-th energy bin of the j-th experiment and 

is the number of bins. Also, let = Yji,j{f^ij ~ ''^ij,HR)^/''^jj,HR5 where 
^i3,'KR. is the predicted average number of events (for AGASA, Fly's Eye and 
Yakutsk) in a model where the HiRes value is the average number of events. 
We find Shi^ = 1.06 with TV = 26 data points, and SModei — 1-20 with N — 2>Q 
data points for all three choices of the absolute energy scale (panels a, b or 
c). The difi^erent experiments are in agreement with each other and with the 
model in the energy range 10^^ eV to lO^'' eV. 

3.2 What is happening above 10^° eV? 

Above 10^° eV, Fig. 2 shows that the Fly's Eye, HiRes and Yakutsk exper- 
iments are in agreement with each other and the model. However, the eight 
AGASA events with energies greater than 10^° eV disagree with the predic- 
tion of the cosmological model (defined by Eq. 1 and Eq. 2), including the 
GZK suppression. The Fly's Eye, Yakutsk and HiRes experiments have a 
combined exposure three times that of the AGASA experiment. The expo- 
sures above 10^*^ eV are, in units of lO'^km^ — yr — sr: AGASA (1.3), Fly's 
Eye (0.9), Yakutsk (0.9), and HiRes (2.2). Together, Fly's Eye, Yakutsk, and 
Hi-Res observe a total of 6 events above 10^° eV (4 events if the Fly's Eye 
energy scale is chosen). 
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Table 1 

Evidence for a GZK suppression. The table compares the number of events expected 
above 10^° eV assuming that there is no GZK suppression with the observed number 
of events for different choices of the absolute energy scale. The expected number 
of events is calculated assuming that the power law J oc £""2.75 dominates 
between 6 x 10^^ eV and 4 x 10^^ eV [25] extends beyond 5 x 10^^ eV. The numbers 
of events are given for the combined exposure of the Fly's Eye, HiRes, and Yakutsk 
experiments. 

Energy Scale Expected Observed 

Fly's Eye 34 4 

AGASA 40 6 

Yakutsk 46 6 

Assuming no GZK suppression, Table 1 compares the expected number of 
events above 10^'^ eV with the number of events observed in the combined Fly's 
Eye, HiRes, and Yakutsk exposure. The differential energy spectrum observed 
by the various experiments at the energy range of 4 x 10^^ eV to 4 x 10^^ eV can 
be fitted by a broken power-law, where the shallower component dominating 
above ~ 6 x 10^^ eV satisfies J oc E-2.75±o.2 ^^^^ ^^^^^ y 43 ^25]). 

The expected number of events in the absence of a GZK suppression was 
calculated by assuming that the cosmic ray spectrum follows the power law 
J oc E~^-'^^ also at energies > 4 x 10^^ eV. Thus there is a > 5(j deficit beyond 
10^° eV relative to the extrapolated lower-energy spectral energy distribution. 
Adopting the steepest allowed slope, J oc E'"^'^^, the expected number of 
events is {21,25,30}, implying a > 3.7a deficit beyond 10^° eV . 



4 Discussion. 

Our most important conclusion is that exotic new physics is not required to 
account for the observed events with energies in excess of 10^° eV, except 
for the AGASA data. Table 1 shows that there is already a strong suggestion, 
> 5(7 (> 3.7(7, depending upon the extrapolated energy spectrum) in the Fly's 
Eye, HiRes, and Yakutsk observations that the expected GZK suppression has 
been observed (see also Fig. 1 and Fig. 2). 

Precision measurements from 10^^ eV to 5 x 10^^ eV are essential for test- 
ing models of UHE cosmic rays, although they are less dramatic than mea- 
surements above 10^*^ eV. At energies > 10^° eV, the predicted number, A^, 
of events in conventional models is uncertain due to the unknown clustering 
scale. To, of the sources, cr (iVpi.edicted)/A^predicted = 0.9(ro/10 Mpc)°-^ [26]. Para- 
doxically, we may need to study carefully cosmic rays with energies below the 
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GZK suppression in order to understand better the origin of the cosmic rays 
beyond the suppression. 
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